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Abstract 

■ 

O |' A Wronskian differential formula, useful for applying the confluent second-order SUSY transfor- 

mations to arbitrary potentials, will be obtained. This expression involves a parametric derivative 
with respect to the factorization energy which, in many cases, is simpler for calculations than the 
previously found integral equation. This alternative mechanism shall be applied to the free particle 
and the single-gap Lame potential. 

a -1 

t^j- ! 1 Introduction 

0^ ' Over time, there has been much interest in the study and generation of exactly solvable potentials in 

quantum mechanics. One of the simplest techniques to carry out these purposes is the supersymmetric 
quantum mechanics (SUSY QM), which is based on the intertwining relationship. In fact, this technique 
together with the factorization method and Darboux transformation are equivalent procedures for gen- 
erating solvable Hamiltonians from a given initial one [TH5]. In the simplest case of this approach both 
Hamiltonians are related by a first-order differential operator. However, the fact that the zeros of the 
transformation function are mapped into singularities of the new potential imposes certain restrictions 
on this method. 

An alternative to avoid this difficulty is to employ the SUSY QM of higher order (for more details and 
applications see [3HH]). The most elementary version of such a generalization is the confluent second- 
order SUSY QM [S1[T3], for which the two involved factorization energies converge to the same value. 
However, the main problem for implementing this method has to do with the difficulty to calculate certain 
integrals which arise in the treatment. 

On the other hand, it has been shown [T4] that in the confluent case the Wronskian formula is preserved 
if solutions closing a Jordan chain of length two are used as seeds for implementing the algorithm. In 
this article we will take advantage of this fact by introducing a differential version of the technique which 
will preserve as well the general Wronskian formula and will avoid to evaluate the previously mentioned 
integrals. In this way, an alternative calculational tool will be available for implementing the confluent 
second-order SUSY QM. 

The paper is organized as follows: in Section 2 we shall review the standard approach to the second- 
order SUSY QM for the confluent and non-confluent cases. A differential version of the Wronskian 
formula for the confluent case will be presented in Section 3. In Section 4 this alternative method will be 



* email: dbermudez@fis. cinvestav. mx 

* email: davidOfis. cinvestav. mx 
t email: jfernandezg@ipn.mx 



1 



Wronskian differential formula for confluent SUSY QM 



Bermudez, Fernandez, Fernandez-Garcfa 



applied to the free particle and the single-gap Lame potentials. A summary of our original results and 
some conclusions are contained in Section 5. 

2 Second-order SUSY QM 

In the second-order SUSY QM one typically starts from the following intertwining relationship [3,5,6 

B+H, (1) 

where 

V(x), (2) 

V(x), (3) 

' (4) 





HB + 




d 2 


H = 




dx 2 




d 2 


H = 




dx 2 




d 2 


B+ = 




dx 2 



dx 



The aim is to determine V(x), g{x) and h(x) supposing that the initial potential V{x) is known. The 
solution to this problem is given by 

V{x) =V+ 2g', (5) 
Kx) = - g - + 9 ^-V + d, (6) 
where g{x) must fulfil the non-linear differential equation 

99" 9' 2 , a ( _, , 9 1 



, j ■ g i [-g' + ^-V + dj+c = o, (7) 

and c, d € K are two integration constants. In order to solve (0 it is used the ansatz [8] 

g\x)=g 2 + 2 19 -2ti, (8) 

where j(x) and £(x) are functions to be determined. Substituting (jSJ) in Equation it turns out that 
^ 2 = c and the following Ricatti equation must be satisfied 

V + 7 2 = ^e, (9) 

with e = d+£. Thus, the initial problem defined by Equations (jT]|4j) has been reduced to find the function 
7 and, consequently, to solve Equation ([9J. This Riccati equation can be linearized by using 7 = u /u, 
which leads to the following stationary Schrodinger equation for H: 

Hu = -u" + Vu = eu. (10) 

The kind of seed solution u employed for constructing the transformation depends on the factorization 
energy e and, consequently, on the sign of c. For c ^ one gets the so called real and complex cases while 
for c = the confluent one is obtained. 

2.1 Non-confluent case (c 7^ 0) 

Let us denote the two different factorization energies by e\ = d + c 1 / 2 and £2 = d — c 1 / 2 , which includes 
both, the real and complex cases. Note that the ansatz ([8]) indeed gives place to two equations 

g' = g 2 + 2 7l .g - ( Cl - e 2 ), (11) 
g' = g 2 + 2 l2 g - (e 2 - e x ). (12) 
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By subtracting both it is obtained 

9 = ^—2- = -{]n[W(u 1 ,u 2 )]y, (13) 
71-72 

where W(f, h) = fh' — f'h is the Wronskian of / and h. The Wronskian in (|13[) should not have zeros 
in order to avoid singularities in g and, consequently, in V. Substituting (113[) in ([5]) one gets 

V = V -2{ln[W( Ul ,u 2 )]}". (14) 

This expression has been used to construct a wide variety of potentials V departing from a given initial 
one V by choosing two appropriate solutions of (fTUl) . In fact, this treatment is so versatile that it can be 
implemented by taking as seeds even the Gamow vectors, which gives place to complex potentials with 
known spectra (see e.g. [15]). 

2.2 Confluent case (c = 0) 

In the confluent case (c = 0) both energies converge to just one, e 2 — > ei, and the ansatz of Equation (jHJ) 
becomes 

g' = g 2 + 2 7 i<7. (15) 
This Bernoulli equation has a general solution given by 

g(x) = -{ln[w(x)}}', (16) 

where x 

w(x)=w - uf(y)dy, (17) 

Jxq 

with wo and xo being real constants which can be chosen at will in order to avoid singularities in g(x). 

On the other hand, let us consider the following pair of generalized eigenfunctions of H, of first and 
second rank, associated to e\ [MlUBlflT] . 

(fl--ei)ui = 0, (18) 
(fl--ei)ua = t*i, (19) 

which is known as Jordan chain of length two. By solving Equation (|19l) for u 2 through the method of 
variation of parameters, supposing that u\ is given, we get 



u 2 = ( k + I -Yr^dx ) ui{x). (20) 



w(x) 
u\{x)' 

Moreover, by using the following Wronskian identity 

W(f,hf) = h'f, (21) 
which is valid for two differentiable arbitrary functions / and h, it is straightforward to show that 

w (x) = W(u u u 2 ). (22) 



Therefore, the Wronskian formula of the non-confluent second-order SUSY QM given by Equation (|14[) 
is preserved for the confluent case [14] . Moreover, it can be used to construct a one-parameter family of 
exactly solvable potentials for each solution u\ of the initial stationary Schrodinger equation associated 
to €\. However, if u\ has an involved explicit form the task of evaluating the corresponding integrals is 
not simple. In the next section we shall present an alternative version of the Wronskian formula for the 
confluent case which will make unnecessary the evaluation of the integrals of Equations ((T7)) and ([20]) . 
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3 Wronskian differential formula for the confluent SUSY QM 

Let us look for now the general solution of Equation (fl9l) in a slightly different way. Let Ux denote 
once again the given solution of (p~8|) . It is well known that the general solution of the inhomogeneous 
second-order differential equation (|19j) takes the form: 

u 2 = u% + v%, (23) 

where u\ is the general solution of the homogeneous equation and denotes a particular solution of the 
inhomogeneous one. Since the homogeneous equation is of second order, it has two linearly independent 
solutions. They can be taken as u\ and its orthogonal function u\ defined by W{u\.u\) = 1. The last 
equation can be immediately solved for , yielding 



dx 
u\{x)' 



ui{x) = Ul (x) I . (24) 



Then, it turns out that 



u£ = Cux + Dui, (25) 
with C, D e M. 

In order to find the particular solution u^, let us suppose from now on that u± and its parametric 
derivative with respect to ex, ^f^, are well defined continuous functions in a neighbourhood of e\. Hence, 
by deriving Equation (|18[) with respect to t\ it is obtained: 

(H-ex)^ = ux, (26) 
oe\ 

where the partial derivatives of u\ with respect to e\ and x have been interchanged. It should be clear 
now that (compare Equations (TT9")) and (|2"B"|) ) 

m 

is the particular solution of the inhomogeneous equation we were looking for. Finally, the general solution 
of equation (JTHJ) is given by 

u 2 = C Ul +Dui + ^. (28) 
oex 

From this equation we can easily calculate the Wronskian of the two solutions of the Jordan chain as 

W(ux,u 2 ) = D + w(ux,^\. (29) 
Thus, the general Wronskian formula of Equation (|14p becomes now 

V = V-2{ln\D + w(ux,^]}} , (30) 



which represents an alternative way to calculate the new potential V through the confluent second-order 
SUSY transformation. 

Note that a special case of equation (|30p has been addressed previously, for D = and the free particle 
potential [T8l[T9] . In these works, the particular solution was taken directly as the seed solution u-i 
and thus the constant D, which arises from the non-trivial term involving the orthogonal function 
(see the second term of the right hand side of Equation (|28p). never appears in those treatments. 

An additional point is worth to remark: without the constant D the confluent second-order SUSY 
partner potential V will often have singularities. The freedom we have here for choosing this constant 
endows us with the possibility to generate families of non-singular potentials for a wide set of factorization 
energies. 
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4 Applications 

We are going to use Equation (|30l) now to implement a confluent second-order SUSY transformation for 
two simple systems. The first of them is the free particle, where both the differential and the integral 
versions of the confluent SUSY QM are easily applicable since the derivatives and the integrals involved 
are not difficult to calculate. The second one is the single-gap Lame potential, for which the previously 
found integral Equation (fTTf is not easy to apply, since the integrals of elliptic functions are complicated 
to evaluate. As far as we know, the confluent second-order SUSY transformation has been never applied 
before to this potential. 



4.1 Free particle 

The free particle is not subject to any force so that the corresponding potential is constant; without loss 
of generality, let us take V(x) = 0. In order to obtain non-singular confluent second-order SUSY partner 
potentials one has to use as transformation function, in general, a solution u\ to the stationary Schrodinger 
equation f| 10[) such that W{u\, U2) ^ V 1 6 R. This is achieved by demanding that u\ vanishes at one 
of the boundaries of the re-domain (see [13l[T4]). In particular, for the free particle these solutions are 
{e KlX , e~ KlX } with the condition that k\ and e\ satisfy the 'dispersion relation' e\ = — ref, «i > 0. 

We are going to use one of these solutions to perform the SUSY transformation, e.g., u\ — e KlX ; the 
other case can be obtained through a spatial reflection. Thus, the parametric derivative can be calculated 
using the chain rule as 

dui d,K\ dui xui xe KlX 
de\ dei dn\ 2k\ 2k\ 

We can evaluate easily the Wronskian of u\ and S^- by using once again equation ([2~l"T) : 

wL,¥)=-£ = -^. (32) 

V dei J 2ki 2ki v ' 

Now, replacing (f3"2"j) into (j3"0)) for calculating the confluent second-order SUSY partner potential V of the 
free particle it is obtained: 

{2Dn 1 -e 2 ^ x Y' { ' 

Due to the 'dispersion relation' (ei = — k\, n\ > 0) there is a natural restriction on the factorization 
energy, namely, t\ < 0. Besides, in order to obtain non-singular transformations the parameter D has to 
be restricted [T31IT4] . Indeed, for u\ = e KlX we have that the non-singular domain is given by D < 0, and 
reparametrizing as D = — (2ni)^ 1 e 2KlX ° , xq € K, we can simplify (j33j) to obtain 

V = -2K?sech 2 [Ki(ir - ar )], (34) 

which is the Poschl- Teller potential with one bound state at the energy Eq = €\ = It is worth 

to note that this result has also been obtained through first-order SUSY QM and by using the integral 
formulation of the confluent case [13] ■ It is plausible that any non-singular SUSY transformation which 
departs from free particle and creates just one bound state leads precisely to a Poschl- Teller potential 
(see also [20]). 

An illustration of a confluent second-order SUSY partner potential V, generated through this formal- 
ism from the free particle, is shown in Fig. [1] 

Note that in some previous works [I8lll9j. the confluent second-order SUSY (or Darboux) transfor- 
mation in this differential version was implemented for the free particle with D — and using another 
transformation function, namely, u\ — sin[fci(a; + xq)} with ei = k\ > 0; however, by doing so one 
will deal only with singular transformations. Following the formalism of this work we have obtained a 
one-parameter family of non-singular potentials for each t\ < 0. 
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Figure 1: Confluent second-order SUSY partner potential of the free particle, obtained through Eq. (|34 
for ei = —1 and xq = 3. 



For the free particle the integral and differential Wronskian formulae have been applied easily, since 
the involved integrals can be simply evaluated. Nevertheless, there are some other potentials for which 
the calculation of the corresponding integrals looks complicated but the differential formalism can be 
applied straightforwardly. We will show next an example of this situation. 



4.2 Single-gap Lame potential 

The Lame periodic potentials are given by [2TM23) : 

V(x) = n(n + l)m sn 2 (a;|m) 



n(n + 1) 



p(x + iK(l - m)) + -(m+ 1) 



(35) 



where sn(x\m) is a Jacobi elliptic function whose real period is T 
elliptic function with 

K (to 



4:K(m) and p(x) is the Weierstrass 



(1 



0)V2' 



(36) 



being half the real period of V(x). The potentials (|35p have 2n + l band edges which define n + 1 allowed 
and n + 1 forbidden bands. They belong to a class of finite-gap periodic systems where the non-linear 
supersymmetry plays an important role. For example, Lame potentials have been used to model a non- 
relativistic electron in periodic electric and magnetic field configurations which produce a ID crystal [24] . 
In addition, these potentials admit isospectral super-extensions |25j and they can be used to display 
hidden symmetries in quantum dynamical problems, specially in soliton dynamics |26j . Note that Lame 
potentials are particular cases of the associated Lame potentials, which have been studied previously in 
the context of higher-order SUSY QM [37]. 

In this work we shall deal with the single- gap Lame potential obtained with n = 1. The spectrum for 
the Hamiltonian associated to this specific potential is given by: 



Sp(H) = [to, 1] U [1 + m, oo), 



(37) 



i.e., it is composed by a finite energy band [to, 1] plus a semi-infinite one [1 + to, oo) (see the white region 
in Fig. [3J. The structure of the resolvent set of H is similar, namely, there is a semi-infinite energy gap 
(—oo, to) plus a finite one (1, 1 + to) (observe the dark zone in Fig. [2]). 

As in the previous case, in order to implement the confluent second-order SUSY transformation we 
are going to use an appropriate seed solution u\ associated to a factorization energy e\ which is inside 
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Figure 2: Spectrum of the Lame potential with n = 1. The white bands correspond to the allowed energy 
region, a semi- infinite [1 + to, oo) and a finite one [to, 1]. The dark region corresponds to the energy gaps, 
a semi- infinite (— oo, m) and a finite one (1,1 + to). 



one of the energy gaps and such that W(u\,U2) ^ V x £ R. For our example this can be achieved by 
choosing u\ as one of the two Bloch functions associated to t\ [2"2"ll23| . i.e., 



o-(a/) cr(a; + <5 + U)') _, 



C(5) 



1 cr((5 + w') cr(a; + cj') 
1//3 _ cr(w') cr(a; — (5 + w') 
1 ~~ cr(-<5 + u/) ct(x + o/) 



(38) 
(39) 



where u> = K(m) and u/ = iK(l — m) are the real and imaginary half-periods of p(x) [28) . and a and £ 
are the non-elliptic Weierstrass functions [29] . 

Note that /3 is defined by the relation u^{x + T) — (3u^(x) and then (3 = exp[2S((uj) — 2uj^(S)}. Besides, 
by expressing it as j3 = e lK , then n = 2i[ll>£(<5) — ^C( w )] ( U P to an additive multiple of 2iri) which is known 
as the quasi- momentum |25j . The displacement 8 and the factorization energy e\ are related by [23 : 

ei = |(m+l)-p(<5). (40) 

In order to calculate (|30[) let us choose the first Bloch function as transformation function, namely, 
u\ = uf. It is worth pointing out that we are using one Bloch state to perform the SUSY transformation, 
even when these states are not normalized. Nevertheless, one of the advantages of the confluent algorithm 
is that it does not require normalized states to perform the transformation. 

We are going to calculate next its parametric derivative with respect to ei, for which we will employ 
the following relationships between cr(x), ((x), and p(x) [25j : 

a'{x) = a(x)C(x), (41) 
C'(x) = -p(x), (42) 



a(2x 

Thus, using the chain rule and Equation (|4"0")) it is obtained 



P'W = (43) 



du\ d5 dui ( dp\ 1 dui 



dd de x 85 \dS J 85' ^ 
An explicit calculation produces: 

^ = {((x + 6 + «,')- C(8 + J) + xp(S)] Ul , (45) 
and thus the Wronskian of Equation (|30)) can be obtained by using once again Equation (l2"Tj) : 
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x 

Figure 3: SUSY partner potential V(x) (top) and the probability density of its associated bound state 
(bottom), generated from the Lame potential Vix) for n = 1. The parameters were taken as m — 0.5, 
ei = 0.1, a; = and D = -45. 



which defines the auxiliary function f(x). 

Finally, from equation (|3T)]) the new potential V can be calculated analytically. 

V _ v , 2 + 4[C(x + 5 + lu') - ((x + J) - OMDu^ + /) 

{Duf + ff ■ 

Two potentials obtained through this method are shown in the top of Fig. [3] and Fig. 2] They 
correspond to two different cases, for which either the factorization energy belongs to the infinite gap or 
to the finite one. Note that the shape of the new potentials (continuous lines) are really different compared 
to the original one (dashed lines), and between them. Indeed, it can be seen that the new potentials are 
in general non-periodic, although they become asymptotically periodic. Note that this periodicity defect 
of V(x) arises due to the creation of a bound state at an energy which coincides precisely with ei (inside 
an initial energy gap). The width and the position of this periodicity defect in general coincides with the 
x-domain where the new bound state 

4?(x)k ^— ^ (48) 

has a non-trivial probability amplitude. For these two cases, the corresponding probability densities 
(x)\ 2 are shown in the bottom of Fig. [3] and Fig. |H 
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Figure 4: SUSY partner potential V(x) (top) and the probability density of its associated bound state 
(bottom), generated from the Lame potential V(x) for n = 1. The parameters were taken as m = 0.1, 
ei = 1.05, x = and D = 20. 

Let us note that a similar physical situation, induced by a non-confluent second-order SUSY transfor- 
mation, was found elsewhere [2"2l2"5] . The main advantage here is that we are using just one seed solution 
to create a bound state inside a given energy gap. Moreover, the explicit expressions obtained from our 
treatment become shorter than those derived by the non-confluent algorithm. Particularly interesting is 
the case in which the factorization energy t\ is inside the finite gap, so that a bound state is created at 
this position. In such a situation, if the non-periodic potential V is perturbed by an additional inter- 
action, the new bound state could be used as an intermediate state to perform transitions between the 
finite energy band and the infinite one. Note that the new bound state of Eq. (|48|) is known as localized 
impurity state in solid state physics (see Ch. 5 of [30]). 

5 Conclusions 

In this article we have introduced a differential version of the confluent second-order SUSY transformation, 
as an alternative to generate new exactly solvable potentials which avoids the need to evaluate some 
integrals arising in the formulation elaborated previously [T31H3] ■ Moreover, we have found a differential 
formula that generalizes the one used in soliton theory |18(I19|. The main advantage rests in the fact 
that families of non-singular potentials can be constructed by appropriately varying the new constant 
D (see Eq. (f30|)). We have successfully applied this technique to the free particle and to the single-gap 
Lame potential. In the last case it was shown that, although the initial potentials are periodic, the SUSY 
generated ones are non-periodic, with a periodicity defect arising due to the creation of a bound state 
inside an initial energy gap. It was proposed that, under certain appropriate circumstances, this bound 
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state could be used as an intermediate state to perform transitions from the lower energy band to the 
infinite one. 
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